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We calculate quasi-normal f- and g-modes of a neutron star with density discontinuity, which may 
appear in a phase transition at extreme high density. We find that discontinuity will reflect largely 
on the f-mode, and that the g-mode could also be important for a less massive star. 
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I. INTRODUCTION 

Gravitational wave interferometers such as LIGO, 
VIRGO, GEO600, and TAMA300 on earth, and LISA in 
space are almost ready to observe directly the gravita- 
tional waves from the Universe. The significance of direct 
observation of gravitational waves is that we obtain new 
method to observe the Universe by gravitational waves 
(the so-called gravitational wave astronomy), and that 
we may be able to verify the theory of gravity, and that 
new physics at high density or high energy region can be 
found by observing a merger of a binary neutron star. 

The one of the sources of the gravitational waves is sup- 
posed to be supernova explosion, which occurs at the last 
phase of a massive star. After explosion, a newly formed 
compact object may oscillate violently, and the emitted 
gravitational waves carry the information of the sources. 
The oscillation will eventually damp out because that the 
gravitational waves carry the energy. The merger of two 
coalescing black holes or neutron stars is another source. 
It is naturally expected that the final object will also os- 
cillate. If a final compact object is a neutron star, its 
structure will reflect on the emitted gravitational waves. 
There should be some relation between a structure of a 
neutron star and the emitted gravitational waves. The 
structure of a neutron star not yet well-known because 
the equation of state is uncertain at the high density over 
the nuclear density. Therefore, although main composi- 
tion of matter in a neutron star is neutron, it is possible 
that there appear new phase in extremely high density 
such as a pion condensation, or a kaon condensation. In 
such a situation, a phase transition may occur, and if 
it is strongly first order, one may expect that there ap- 
pears the discontinuity of density at some density as the 
case of the kaon condensation ||]. If we know in advance 
the relation between the structure of a neutron star and 
the emitted gravitational waves, we will get the informa- 
tion or the restriction on the equation of state over the 
nuclear density. In particular, if we find any specific fea- 
ture in the emitted gravitational waves from a neutron 
star with discontinuity, we may be able to predict the 
existence of a phase transition in a neutron star. This is 
the main subject of the present paper. We will just focus 
into quasi-normal modes of a neutron star. 

Although the gravitational waves emitted in a super- 



nova explosion or from a merger of a binary may show 
very complicated wave form, only quasi-normal modes 
will eventually remain. We expect that these quasi nor- 
mal modes contain a lot of information about a neutron 
star. Hence, it is very important for us to study those 
modes for a possible situation. The quasi normal modes 
are classified into various types. For the modes coupled 
to the fluid of a non-rotating neutron star, we have f-, p- 
and g-modes. 

The damping rate of these modes is very small com- 
pared with the oscillation period of fluid because mat- 
ter couples to gravity very weakly. There exists only 
one f-mode for each I. The p-mode which is caused by 
the pressure of fluid. The g-modc is originated by den- 
sity discontinuity or temperature of a star. Its damping 
rate is usually quite small, so few calculation has been 
done. There also exists another mode (the so-called w- 
mode) which is related to metric perturbations of space- 
time rather than fluid oscillation. This mode has high 
damping rate in comparison with modes related to fluid 
oscillation. 

For the case without density discontinuity, the rela- 
tion between quasi normal modes of a star and equation 
of state has been studied by many authors |3|, Q, [||. 
For g-mode, there are a few works. Newtonian case is 
studied in |^ and also in which includes the effect 
of temperature of neutron stars by use of the Cowling ap- 
proximation 1^. For the relativistic case, Finn adopted 
an approximate method which is called the "slow motion 
formalism" and calculated the imaginary part by the en- 
ergy loss for the case of density discontinuity 0, 

Although these works are very important, they have 
not discussed density discontinuity at high density, which 
may appear via a phase transition. In this paper, we shall 
study a neutron star with such a discontinuity, and in- 
vestigate its quasi normal modes. Just for simplicity, we 
assume that equation of state is given by a polytrope 
(P = i^Tp-^+^Z"). We discuss only even parity modes, 
which are coupled to oscillation of matter fluid in a neu- 
tron star. We then assume that there exists density dis- 
continuity due to a phase transition and density discon- 
tinuity. 

First we present the basic equations to solve in § p| in- 
cluding the equation for a background neutron star with 
density discontinuity. In §111, we show an explicit struc- 
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ture of a neutron star with the density discontinuity and 
numerical results for quasi normal modes. We also com- 
pare our results with those obtained by the relativistic 
Cowling approximation [pl| . Conclusion is given in §IV. 
In Appendix ^ we present the explicit boundary condi- 
tion at the center of a star. In Appendix we briefly 
summarize the method of how to calculate the quasi nor- 
mal modes. 

In this paper we adopt the unit of c = G = 1, 
where c and G denote the speed of light and the gravita- 
tional constant, respectively, and the metric signature of 

(-,+,+,+)• 



II. BASIC EQUATION 



As for the equation of state, just for simplicity, we 
adopt a polytropic equation of state except at a discon- 
tinuous surface of density, which is assumed to exist in- 
side of the star. Then, the equation of state is given by 
two equations; 



for < r < i?dis, (2.7) 
for i?dis <r<R, (2.8) 



where i?dis is the radius of a discontinuous surface. The 
values of two adiabatic indecies (f^(+) and '^(-)), and 
those of two coefficients (-ft^(+) and may be dif- 

ferent. However, since the pressure should be continuos 
at the discontinuous surface i?dis, we have one constraint 
on those parameters. For example, is fixed as 



A. stellar model with discontinuity 

We first construct a neutron star model with disconti- 
nuity of density to calculate its quasi-normal modes. A 
static and spherically symmetric spacetime is described 

by 
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92 + sin2 



where $,A are metric functions with respect to r. 
mass function m{r), which is defined as 



(2.1) 
A 



m(r) 



1 



(1 



-2A 



satisfies 
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The equilibrium condition of a stellar model is given by 
the Tolman-Oppenheimer-Volkoff equation. 



dP _ ip + P) (m + Attt^P) 
dr r (r — 2m) 



(2.4) 



where the function p and P are the energy density and 
the pressure of the fluid. The potential $ is given by 



dr 



m + 47rr^P) 
r (r — 2m) 



(2.5) 



With an additional equation, i.e. the equation of state, 
we can solve these equations. 

As for the boundary conditions, there are two bound- 
aries; the center of a star and the surface. The stellar 
radius R is determined by the condition that the pressure 
vanishes at the surface, i.e. P{R) = 0. The total mass 
M is given by m{R), and the potential $(r) is smoothly 
connected at the surface, i.e. 



= 2<E'(fl) 



1 



2M 



(2.6) 



[Pdis) 



(2.9) 



Now we have six parameters; the central density of the 
star the inner density at the discontinuous surface 
pI^s ' ratio of the outer density to the inner density 

at the discontinuous surface Pdis' /Pdis' > the inner poly- 
tropic index "■(+), the outer polytropic index and 
the outer polytropic coefficient K(^_y Some examples 
of this equation of state with discontinuous density are 
given in Figure |l|. 
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FIG. 1. Equation of state with density discontinuity. 
The density and pressure are normalized by the stan- 
dard nuclear density pnuc ~ 2.68 x 10^* (g/cm^). We 
set pc = 3.0 X 10^^ (g/cm^ 
pi7sVpit'=0.7, 7^(_) = 100(km^ 



Pdis 



1.2 X lO^'^ (g/cm=*) 
and n(_) — 1.0. 



At the center, a regularity condition should be imposed. 
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B. perturbation equations inside the star 



He re w e present the basic equations for perturbations against the background metric of a spherically symmetric 
star (2.1). The perturbation is given by 



(2.10) 



The perturbations against a spherical background can be decomposed by the spherical harmonics Yj^, and are classified 
into two by its parity, i.e. odd or even. Because we are interested in to discuss relation between a discontinuity in 
stellar density and the gravitational waves emitted from such a neutron star, in this paper, we only study even-parity 
perturbations, in which gravitational waves are coupled to fluid oscillation. 
Following the method by Lindblom and Detweiler , hf^^, is described as 



i TT„2<S> 
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iujr^+^Hi r^He^^ 

r^+^K 
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I r^He 



(2.11) 



where H,Hi, and K are the perturbed metric function with respect to r. The fluid perturbations are described by 
the Lagrangian displacement vectors ^* = (^'", ■C^, C"^) as follows: 
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where W and V arc functions with respect to r. Furthermore we introduce the new variable X instead of V to make 
the boundary condition at the surface of a star simple (see below); 
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The perturbation equations of the Einstein equations are given by 



dHi 
dr 



dK 

dr 

dW 
dr 

dr 



; + 1 + e^^ 

r 



47rr^(P - p)e 



2A 



,2A 



r 

/(^ + 1) 
2r 



H + K + 1Qt:{P + p)V 
1 



Hi 



-H - 



r dr 



K 



{P + p)e^W, 



-W + re' 



1 



X 



2 



--X + {P + p)e'' 
r 



1 1 



1 



- ]H-- \ uj^re 



2r 2 dr J 



dr r 
l{l + l) d<S> 

j,2 flj. 



.2„„-2$ 



^(/ + 1) 

2r 



V 



An{P + p)e^-r-{^[-e'^'^-f] 
dr \r^ dr J 



Hi 



W 



3m l{l + l) 



Airr-^P 



+r'e-^^ 



■ 2 -2* Kl + 1) d<S> 



H - Snr^e^'^X 
Hi 



2r dr 



(2.16) 
(2.17) 
(2.18) 



(2.19) 



3 
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where 7 is the adiabatic index of the unperturbed steUar model defined by 

p + P fdP" 
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Eqs. (Pl)~(|2l9|) give a set of differential equations for the variables Hi, K^W and X, while Eqs. ( |2.20| ) and ( |2.2l| ) 
are the algebraic equations for the variables H and V . 



C. the perturbations outside the stars 

In the region outside the star, the perturbations are described by the Zerilli equations. Setting m — M, W — V — 0, 
and replacing Hi,K with new variables Z, dZ/dr^ defined by 



, " A(A + ly + 3XMr + 6Af ^ ^ dZ 

r K — — Z H 

r2(Ar + 3M) dr,' 

, , , . Ar2 - 3AAfr - iM^ r^ d.7. 

r^+^Hi = r^^rr. 



(r - 2M) ( Ar + 3M) r - 2M dr, ' 
we recover the Zerilli equations. Here r* is the tortoise coordinates and A = Z(^ + l)/2 — 1. 



(2.23) 
(2.24) 



D. boundary conditions 



In order to solve the above perturbation equations, we have to impose the boundary conditions. At the center of 
the star, we have regularity condition for the perturbative variables. We expand all variables by Taylor power-series 
near r = in a similar way of Detweiler and Lindblom (see Appendix ^ . Another boundary condition is that 
the pressure vanishes at the stellar surface. Then the perturbed pressure X, AP is given by 



'X 



(2.25) 



also vanishes there. This boundary condition at the surface of star is equivalent to X{R) = 0. The last condition 
is that there exist only outgoing gravitational waves at infinity. We also have a junction condition at the surface 
of discontinuity, where the variables of Hi,K,W, AP (and then X) must be continuous. Note that, H and V are 
determined by Eqs. (^.20) and (2.21). H is continuos but V is discontinuous. 



III. NUMERICAL RESULTS 
A. models of neutron star 

As a concrete example of density discontinuity, we know kaon condensation which may occur at p^^J = 1.466 x 

10^^ g/cm'^ with the discontinuity PjisVPdis^ — 0.414 for a cold neutron star. Hence we consider density discontinuity 

around this value, i.e. we study the following three cases: p^^J = 8.0 x 10^^, 1.2 x 10^^ and 1.6 x lO^^g/cm'^. First 
we briefly discuss a models of a neutron star with density discontinuity. As mentioned before, we have six unknown 

parameters: 't-(-), -f'^(-), '^(+), Pdt5\ PdisVPdts^ ^^^d pc- We set — 1.0 and = 100 (km^), which are usually 

adopted in a neutron star model. The rest four are «(+), Pd^\ PdisVPdis^ pc- 

First we show the gravitational mass M with respect to a central density pc for several values of Pdi^\ PdisVPdis^ 
in Figs. 
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FIG. 2. The relation between the central density and the gravitational mass of a neutron star M for n(_|.) = 1.0. Fig. 
2(a) is the case of p^ts' = 8-0 ^ 10^'' (s/cm^) , while Fig. 2(b) is the case of p^+' = 1.6 x 10^^ (g/cm^). The numbers associated 
with each curve are the values of Pjis' /PdTs' ■ 




As for inner polytropic index, we set n(+) = 1.0 (Fig. ||), 0.7 (Fig. ||), 1.5 (Fig. ||). For each polytropic index, we 
consider two models with different density at discontinuity ip^t,-* = 8.0 x 10^"* (a) and 1.6 x lO^^g/cm'^ (b). We depict 
the M-pc diagram for each value of 9"^^ I P'asI ^ ^-^ (^o discontinuity), 0.9, 0.7, 0.5, and 0.3. Fig. || for — 1.0 
shows that the maximum masses for each value of 1 Pa\s '^'^^ depend on density discontinuity p^t-' . As 1 P]x2 
gets small, the maximum mass decreases and the stable region {dM / dpc > 0) becomes narrower and moves to high 
density region. 

For «(+) 7^ 1.0, the maximum mass depends on p''^^ ■ In fact, for «(+) < 1.0 (Fig. the maximum mass decreases 

as p^^' increases, while for rt(+) > 1.0 (Fig. the maximum mass increases. As the equation of state of inner region 
becomes more stiff, i.e. as 7i(_|_) decreases, the maximum mass gets large. The other properties about maximum mass 



5 



for n(+) = 0.7 and 1.5 are same as those of n(+) = 1.0. Note that for the case of = 1.5, all stellar models are 
unstable if discontinuity gets large. 

In what follows, we discuss only stable stellar models, i.e. dM/dpc > 0. In our analysis, we fix the mass of a neutron 
star AI (1.2M0, 0.5Mq) because M may be determined by other observations. In Table |, Table |l| and Table 
summarize the properties of background stellar models. 



Ill 



we 




B. Quasi- normal modes for a stellar model without discontinuity 



In order to find a specific property in quasi-normal modes of a neutron star with discontinuity, first we calculate of 
I = 2 in the case of a stellar model without discontinuity as reference. In this case, we have to set two parameters; the 
polytropic index n and polytropic coefficient K, if we fix the total mass M. We adopt the following values; (n, K) = 
(1.0,100km2), (1.0,200km2), (1.5, ISkm^/^), (2.0,3.5km) for M = 1.2Mq. We also the case of M = 0.5Mq with 
(n, K) = (1.0, lOOkm^) to examine the mass dependence. The properties of the neutron star models are summarized 
in Table 

The quasi-normal modes of the gravitational wave emitted from such stars are plotted in Figure |^ and listed in 
Table ^. The plot point corresponds to the so-called f-mode of each stellar model. In this table we show the mode 
frequency u by two normalizations: one is wM and the other is u){R^/My/^. Although the f-mode is determined by 
the radius as well as the mass, we may not know the radius from observation. Then we can plot the data by uiM. In 
the figure, however, we only show uj{R^ /M)^^'^ . In a neutron star, we usually find many p-modes as well as f-mode, 
but we only show the lowest f-mode here because that we are interested in new mode which appears in low frequency 
for a star with discontinuity. 
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FIG. 5. The quasi normal mode (f-mode) of a neutron star without density discontinuity. The labels associated with each 
mark correspond to the models in Table 
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From Fig. |5[ we find that the real part of eigenfrequency gets large as equation of state is softer. For example, the 
frequency in the case 'd4' is by 47% larger than that in the case 'dl'. However, this tendency becomes reverse if we 
use another normalization (wM) because the radius of a neutron star constructed by soft matter is larger than one 
constructed by stiff matter (see Table IV). We also find that the imaginary part of eigenfrequency increases as the 
stellar mass increases, but the real part is insensitive to the mass. 



C. Quasi- normal modes for a stellar model with discontinuity 



In the case with discontinuity, we have three unknown parameters: «(+), p)^2 ^^"^ Pa\1 I P^a\1 ■ First we show the 
eigenfrequency of quasi-normal modes of Z = 2 for the case of nj^) = 1.0. We fix the mass of a neutron star as 

M = I.2M0 or O.SMq. Now two free parameters are p"^2 ^'^^ P'asI I P">x2 ■ 

We show the f-mode for M = I.2M0 and O.SM© in Fig. ^. The numerical values are given in Table VI and Table 



VII. 
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FIG. 6. The quasi normal mode (f-mode) of a neutron star with density discontinuity. We set nj^.) = 1.0. The left figure 
(a) shows the result of a star with M = 1.2Mq, while the right one (b) corresponds to a star with Af — 0.5Mq. The labels 
associated with each mark correspond to the models in Table | and in Table ^ The data of 'dl' and 'd5', which are the models 
without discontinuity for n = 1.0 and K — 100 (km^), are also plotted as a reference. 



For f-mode, we find that the mode frequency for less massive star (0.5 Mq) is more sensitive to discontinuity than 
that for massive star (1.2 Mq). For example, the difference of the model 'b5' from the model without discontinuity is 
larger than 28% for 0.5Mq. It is also found that the imaginary part of eigenfrequencies in the model with discontinuity 
is larger than one without discontinuity, and this difference becomes larger as M/R increa ses. 

The results of g-mode for M = 1.2Mq and O.5Af0 are given in Fig. I and in Table ^m] and Table For g-mode, 
we also find the similar properties to the f-mode. The mode frequency in a less massive star is more sensitive than 
that in a massive star, in particular the imaginary part changes by the factor 10^ — 10^ for the same p^^J depending 
on the amount of discontinuity (see Fig.^. Furthermore, the existence of g-mode is due to discontinuity, and it turns 
out that its imaginary part of g-mode is not always so small if the discontinuity is very large. We conclude that the 
damping rate of g-mode may not be ignored in a particular situation. 

Next we calculate the eigenfrequencies of f-mode and g-mode for — 0.7. The results of f-modes and g-mode 
are shown in Fig. |8|(a) with the case of the model without discontinuity, and in Fig. ^(b), respectively. The numerical 
data are given in Table ^ and Table XI, respectively. By comparing the upper panel in Fig. ^(a), we find that the 
effect of change of equation of state is quite small. On the other hand, the imaginary part of g-mode becomes by 10 
times larger than the case of n = 1 due to the existence of the inner stiff matter. Furthermore, in this case, even if 
the density is continuous (PdisVPdis^ = l-O); one may expect the appearance of g-mode because of discontinuity of the 
equation of state (?i(+) = 0.7 and = 1.0). However, we could not find g-mode in our calculation. This may be 
because the g-mode does not exist or is too small to be found. 
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FIG. 7. The quasi normal mode (g-mode) of a neutron star with density discontinuity. The setting is the same as Fig. 6. 
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FIG. 8. The quasi normal modes of a neutron star with density discontinuity ((a) f-mode and (b) g-mode). We set n(_|.) = 0.7 
and M = 1.2Mq. The labels associated with each mark correspond to the models in Table |l|. The data of 'dl', which is the 
model without discontinuity for n — 1.0 and K — 100 (km^), are also plotted as a reference. 



D. Comparison with the Cowling approximation 

Finally, by using relativistic Cowling approximation [ ]lT| , we shall calculate the f- and g-modes, and compare 
them to our results. In the Cowling approximation, metric perturbation is ignored and only stell ar os cillation is 



taken into account. Hence, the eigenfrequency has only the real part. The result is given in Table XII for — 
8.0 X 10^'* (g/cm^). The error (the difference from our numerical value) estimated by 

ujcow ~ Tleiuj) 

error = , (3.1 

Ke{uj) 

is also shown. Here wcow is the eigenfrequency f ound by the Cowling approximation, and (^^ow ^^'^ '-^Cow denote 



those of g- and f- modes, respectively. From Table XII, we find that the Cowling approximation is very accurate for 



g-mode if is at least in this region, although it is not so good for the f-mode. 



IV. CONCLUSION 



In this paper we calculate the eigenfrequencies of f-mode and g-mode to examine the dependence of the discontinuity 
of the density if there exists. We investigate with the various discontinuity, by changing p[^J and /OdLVPdk\ foi" the 
case that the total mass is fixed 1.2Mq for n(_|_) = 1.0, 0.7 and O.5M0 for = 1.0. In consequence we find that the 
dependence of the eigenfrequencies of f-mode on the discontinuity of the density is more strongly for the case that the 
total mass is smaller. The effect that the inner matter is stiff, n(^_^_■^ ~ 0.7 and = 1.0, is quite small for f-mode. 
And it is found that the imaginary part of the eigenfrequencies of g-mode become a little large for the case that inner 
matter is stifi^. In general imaginary part of g-mode is quite small, but we find that in particular situation such that 
an inner polytropic index is the same as an outer one and total mass is rather small, the imaginary part of g-mode 
may not be ignored. 

For the considering the circumstances mentioned above, it might be possible to judge whether the discontinuity of 
the density exist or not by the direct observation of the gravitational waves, which might be practicable in future 
observation, for even f-mode in the particular situation which the total mass is small. Of course, if the gravitational 
waves of g-mode is directly observed for such stars as the small massive ones, it is found that the discontinuity of the 
density exist. So we can judge the existence of the phase transition of the density in the neutron stars by the direct 
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observation of the gravitational waves, though we need to collect the detailed data of the eigenfrequencies for some 
characteristic mode. 

Finally we find also that Cowling approximation is quite accurate for g-mode. 
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APPENDIX A: THE TAYLOR EXPANSION AT THE CENTER OF A NEUTRON STAR 



Here we present the explicit perturbation equations (2.16)— (2.21) by the expansion at the center of a neutron star. 
p, P and <& are expanded at r = as 



p(r)=p(") + ip(2)r2+0(r4), 

P{r) - + ipf^'r^ + 1^(4)^4 ^ o(r6), 
\ J 2 4 ^ " 

$(r) = $W + i$(2)^2 ^ 1^(4)^4 ^ o(^6)_ 

\ / 2 4 



where the coefficients are obtained from Eqs. (2_^) — (2^) and ( 2.22 ) as 



<I,(2)^l!Lp(0)+4^p(0)^ 

3 

p(2)= _ (^P(O) +p(0)^ $(2)^ 



p(o)^(o) 



$(4)^ ^^(2) + 2^P(2) ^ 32^^(0)2 ^ 32712 ^(o)p(O)^ 

5 9 3 

p(4)^ _ l^p(O) _^ ^(0)^ ^(4) _ 1 |^p(2) _^ ^(2)^ ^(2)_ 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 

(A6) 

(A7) 
(A8) 



7^°) is the adiabatic index at r = 0, which is defined by the equation of state. The perturbation variables Hi, K, W 
and X are expanded at r = follows, 



Wir)= + iw^(2)^2 ^ (^(^4)^ 

l(r)=X(") + ix(2)r2 + 0(r4), 



1 



r(2)^2 



(A9) 
(AlO) 
(All) 
(A12) 



where the leading order coefficients in these equations are fixed by the perturbation equations ( 2.16| ) — (2.21) as 



1k(0) + f l![^(0) +4^p(0) „ ^^^^-2*'0, ^ ^(0) 

2 \ 3 / 



(A13) 
(A14) 



In the similar way, by using the equations (2.16)— (2.21), we find the second order coefficients in the expansion as 



2 L it I 1 ) 



1^(2; + 3) - 4^P(") I - ^(P(2) + p(2))W^(") + 87r(p(°) + p(°))Qi + Iqo, 



Z + 2 



^(2) _ i([+i)^^(2) _ 4^^p(0) ^ ^(0))^(2) 



2 ' 4 
47r 



p(2)+^(2)+^^(0)(p(0)^^(0)) 



^(0)+(^|[^(0)+4^p(0)j ^(0)+lg^^ 



3*<°'(p(0)+p(0)) 



i^e-2*<'"+27r(p(°)+p("))-^ 



(2) 



(A15) 



(A16) 
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(2) 



2 



+ i±^X(2) + e*™ (P(0) + p(0))jy, 

2 8 



+ 



2 - 2 

- 27r(P(2) + p(2)) _ l^p(0)(p(0) + ^(0)) _ l!L^(0)^(2) 

3 3 

+ 2V +P(0)+p(0)y)^ ' 

(0) , JO)-,] Ti/(2) 



+ $(4) + ^2g-2*(0) j^^(2) _ |I^(0)^ 1^(0)^ ^ ^ 

1 rp(2) + ^±i^^2g-2*""(p(0)+^(0))" 



^(2) _ lg-*(0)^(2) 

2 



_1(P(0) +^(0))^(2) _ 1 Jp(2) ^ ^^c.-e - ^/^-^ vr-' - -e - ^ 

= _le-*'"'$(2)x(°) + i(p(2) +p(2))i^(o) + l(p(o) + p(o))Q^ 

+ [p(4) _ ^p{0)p{2) ^ ^^-2^(0, 1 l^p(2) ^^(2)) _ (p(0) +^(0))^(2)| 

3/12 J 

1 .2^-2*<o)^p(0)+^(0))g^^ 



where Qo and Qi are defined by 



3 / 

:-2*'°)_;(; + i)(^|[^(0)+2,rPw)| 
^(O)p(o)^ ^ 2 3^. 



(0)' 



H and are also expanded at r = as 

1 



t/(r) 



= y(0)+ ly(2)^2^0(r4)^ 



where the 



coefficients are given by 

l/(o)=_lw^(o)^ 

^(2)=Qo + i^(2), 



(A17) 



(A18) 



(A19) 
(A20) 

(A21) 
(A22) 

(A23) 
(A24) 
(A25) 
(A26) 
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APPENDIX B: HOW TO CALCULATE THE 
QUASI NORMAL MODES 



To calculate the quasi normal modes of a neutron star, 
we follow the method developed by Leins, Nollert and 
Soffel ||T3]. First, we change the variables from Zerilli's 
one (Z) to the Regge- Wheeler's one (Q), because the 
effective potential in the Regge- Wheeler equation is very 
simple. The relation between those two variables given 

by B 



{k + 2iujP) {k + 2/3^ f{r)) Q + 2(3 
{k - 2iLuP) g = (k + 2P^f{r)) Z ^2(3 
where 



dQ_ 

dZ 

dr^, ' 



fir 



(3 = 6M, 
K = 4A(A + 1), 
r~2M 



2r^{Xr + 3M)' 



(Bl) 
(B2) 



(B3) 
(B4) 

(B5) 



To find the solution of the Regge- Wheeler with appropri- 
ate boundary condition at infinity, we expand the vari- 
able Q at a finite radius r = as 



Qir] 



\2M 



n=0 



(B6) 



By substituting this into the Regge- Wheeler equation, we 
obtain a four-term recurrence relation for a„ (n. > 2) ; 



+ Pnan + InO-n-l + '5na«-2 = 0, 



where 



2M 



ra ) 
2 I iuira + 



n{n + 1), 



(B7) 

(B8) 
(B9) 



6A/ 



7n = 1 i I n(n - 1) + l{l + 1), (BIO) 

Ta J Ta 

5^ = —{n-i){n+l). (Bll) 

The coefficients a_i, ao, and ai are given by the value of 
Q and dQ/dr ai r — ra as 



a„i = 0, 

Qira) 



ao 



ai 



X{ra) 



X{ra) 



dQ{ra) 
dr 



+ Z ^Q(^a) 



2M 



(B12) 
(B13) 

(B14) 



where 



^^^^ = (2F - 1 



-2iMuj 



(B15) 



The convergence condition (|a„+i/a„| < 1) is now ra > 
4A/ If i? > 4A/, we set = R, while if i? < 4Af, 
we integrate the Regge- Wheeler equation from r = R to 
r = ra{> 4:M) to find the values of Q{ra) and dQ{ra)/dr. 
Furthermore, defining new coefficients d„, /?„, and 7„ as 



di = ai, /3i = /3i, 71 = 71, 
and for n > 2 



7n = In 



7n-l 
ln-1 



(B16) 

(B17) 
(B18) 

(B19) 
(B20) 



the above four-term recurrence relation is reduced into a 
three-term recurrence relation 



d„a„+i -I- /3„a„ -I- 7„a„_i = 0. 



(B21) 



Using this three-term recurrence relation, we find the re- 
lation between Q!„, /9„ and 7„ as the form of continued 
fraction: 



oi _ -7i ai72 ^273 
ao /3i- /32- /33- 

which can be rewritten as 

do7i di72 d273 



= /3o 



/3i- 132- p3- 



(B22) 



(B23) 



where /3o = ai/ao, ctQ = ^1- The eigenfrequency cj of a 
quasi normal mode is obtained by solving the equation 
/(u;) = 0. 



[1] K. S. Thorne, in Black Holes and Relativistic Stars, 
edited by R. M. Wald (University of Chicago Press, 
Chicago, 1998), p. 41. 

[2] T. Tatsumi, and M. Yasuhira, Nucl. Phys. A 670, 218 
(2000). 

[3] N. Andersson, and K. D. Kokkotas, Phys. Rev. Lett. 77, 
4134 (1996). 

[4] N. Andersson, and K. D. Kokkotas, Mon. Not. R. Astron. 

Soc. 299, 1059 (1998). 
[5] K. D. K okkotas, T. A. Apostolatos, and N. Andersson, 

preprint |gr-pc/990107^ . 
[6] P. N. McDermott, Mon. Not. R. Astron. Soc. 245, 508 

(1990). 



13 



[7] L. S. Finn, Mon. Not. R. Astron. Soc. 222, 393 (1986). 
[8] L. S. Finn, Mon. Not. R. Astron. Soc. 227, 265 (1987). 
[9] T. G. Cowling, Mon. Not. R. Astron. Soc. 101, 367 
(1941). 

[10] T. E. Strohmayer, Astrophys. J. 417, 273 (1993). 

[11] P. N. McDermott, H. M. Van Horn and J. F. Scholl As- 
trophys. J. 268, 837 (1983). 

[12] L. Lindblom and S. Detweiler, Astrophys. J. Suppl. Ser. 
53, 73 (1983); S. Detweiler and L. Lindblom, Astrophys. 
J. 292, 12 (1985). 

[13] M. Leins, H. -P. NoUert, and M. H. Soffel, Phys.Rev.D 
48, 3467 (1993). 

[14] S. Chandrasekhar, in The Mathematical Theory of Black 
Holes (ClarendonPress, Oxford, 1983), chap. 4. 



[15] O. Benhar, E. Berti, and V. Ferrari, preprint ^r- 



qc/9901037 . 



14 



TABLE I. The properties of a neutron star constructed by the equation of state with discontinuity. We set M = 1.2Mq and 
n(+) = 1.0. 



label 


piV (g/cm^) 




Pc (g/cm^^) 


R (km) 


Rdis/R (%) 


Mdi./M (%) 


M/R 


al 


8.00 X 10"* 


0.9 


5.017 X 10'" 


7.577 


81.42 


89.23 


0.2339 


a2 


1.20 X 10^^ 


0.9 


4.716 X lO''^ 


7.800 


73.13 


78.20 


0.2272 


a3 


1.60 X 10^^ 


0.9 


4.278 X 10^^ 


8.112 


64.02 


62.90 


0.2184 



TABLE II. The properties of a neutron star constructed by the equation of state with discontinuity. We set M = O.SMq 
and n(+) = 1.0. 



label 


pits' (g/cm^^) 




pc (g/cm^') 


R (km) 


RmJR (%) 


M^iJM (%) 


M/R 


bl 


8.00 X 10'" 


0.4 


1.759 X 10'" 


3.925 


86.99 


98.76 


0.1881 


b2 


8.00 X 10" 


0.5 


5.800 X 10^-' 


6.050 


73.65 


91.10 


0.1220 


b3 


8.00 X 10" 


0.6 


2.484 X lO''^ 


8.465 


56.76 


65.25 


0.08721 


b4 


1.20 X 10^^ 


0.4 


1.733 X 10^® 


4.137 


80.71 


97.16 


0.1785 


b5 


1.20 X lO^'^ 


0.5 


5.127 X lO^'^ 


6.985 


59.38 


76.64 


0.1057 


b6 


1.60 X lO''' 


0.4 


1.680 X 10"^ 


4.407 


74.28 


94.69 


0.1675 



TABLE III. The properties of a neutron star constructed by the equation of state with discontinuity. We set M = 1.2Mq 
and n(_|_) = 0.7. 



label 


pit' (g/cm=^) 


P'.-J/P^' 


Pc (g/cm^) 


R (km) 


Rdi./R (%) 


Mdis/M (%) 


M/R 


cl 


8.00 X 10" 


0.7 


4.237 X 10'" 


7.351 


85.64 


93.61 


0.2410 


c2 


8.00 X 10" 


0.8 


2.787 X 10'^ 


8.488 


79.57 


85.53 


0.2088 


c3 


8.00 X 10" 


0.9 


2.029 X lO"' 


9.387 


73.47 


74.03 


0.1888 


c4 


8.00 X 10" 


1.0 


1.580 X 10'^ 


10.04 


67.38 


60.31 


0.1765 


c5 


1.20 X 10^^ 


0.7 


5.924 X 10'^ 


6.724 


83.11 


91.86 


0.2635 


c6 


1.20 X 10^^ 


0.8 


3.483 X 10'^ 


8.090 


73.88 


78.54 


0.2190 


c7 


1.20 X 10^"^ 


0.9 


2.449 X lO'-""' 


9.046 


65.08 


60.68 


0.1959 


c8 


1.20 X 10'^ 


1.0 


1.903 X 10'" 


9.634 


56.62 


42.14 


0.1839 


c9 


1.60 X 10'' 


0.8 


4.068 X 10'^ 


7.862 


68.95 


71.75 


0.2254 


clO 


1.60 X 10"^ 


0.9 


2.740 X 10"^ 


8.912 


56.79 


46.69 


0.1988 


cll 


1.60 X 10'^ 


1.0 


2.146 X 10"^ 


9.421 


45.71 


25.46 


0.1881 



TABLE IV. The properties of a neutron star constructed by the equation of state without discontinuity. The total mass M 
is fixed I.2M0 or Q.5Mq. 



label 


Mass (Mq) 


n 


K 


Pc (g/cm-^) 


R (km) 


M/R 


dl 


1.2 


1.0 


100 (km^) 


2.427 X 10'" 


9.296 


0.1906 


d2 


1.2 


1.0 


200 (km^) 


5.507 X 10" 


15.10 


0.1173 


d3 


1.2 


1.5 


15 (km*/3) 


2.691 X 10" 


23.49 


0.07542 


d4 


1.2 


2.0 


3.5 (km) 


1.153 X lO'"* 


38.79 


0.04568 


d5 


0.5 


1.0 


100 (km^) 


5.067 X 10" 


11.57 


0.06384 
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TABLE V. The numerical data of quasi normal mode (f-mode) of the models in Table IV 



These data are plotted in 



label 


7^e (ujM) 


Im {loM) 


7^e 1 




Tm {u^jR'^/M^ 


dl 
d2 
d3 
d4 

d5 


0.09784 
0.04869 
0.02990 
0.01688 
0.01974 


4.643 X lO"'' 
1.145 X 10-^ 
3.572 X 10-*' 
8.075 X lO-'^ 
1.374 X lO"** 


1.176 
1.211 
1.443 
1.729 
1.224 


5.579 X 10"'' 
2.849 X 10-" 
1.725 X 10"'' 
8.270 X 10-^ 
8.517 X 10"^ 


TABLE VI. The numerical data of quasi normal mode (f-mode) of the models in Table] 
These data are plotted in Fig.^ (a). 


] (M = I.2M0 and n(+) = 1.0). 


label 


Tie (ljM) 


Im {ujM) 


Tie 1 




Tm {u^R^/M^ 


al 
a2 
aS 


0.1317 
0.1273 
0.1209 


6.856 X 10-" 
6.604 X lO-"^ 
6.221 X lO-"^ 


1.164 
1.176 
1.184 


6.062 X 10-" 
6.098 X 10-'' 
6.093 X 10-'' 


TABLE VII. The numerical data of quasi normal mode (f-mode) of the models in Table|^ (M = O.SM© and = 1.0). 
These data are plotted in Fig.^ (b). 


label 


TZe (ljM) 


Tm {ujM) 


TZe 1 




Tm {u^R^/M^ 


bl 
b2 
b3 
b4 
b5 
b6 


0.1078 
0.06266 
0.03767 
0.1060 
0.05416 
0.1025 


5.384 X 10"' 
1.916 X 10-^ 
5.867 X 10-** 
5.185 X 10"^ 
1.253 X 10"^ 
4.780 X 10-^ 


1.321 
1.470 
1.462 
1.406 
1.576 
1.495 


6.600 X 10-" 
4.495 X 10-" 
2.278 X 10-" 
6.878 X 10-" 
3.645 X 10-" 
6.971 X 10-" 


TABLE VIII. The numerical data of quasi normal mode (g-mode) of the models in Table^ (M = I.2M0 and n(+) = 1.0). 
These data are plotted in Fig.|^ (a). 


label 




Im {ujM) 


TZe 




Tm (ujy/R^/M^ 


al 
a2 
aS 


3.090 X 10"^ 
3.598 X 10-2 
3.884 X 10-2 


4.410 X 10-'^ 
1.955 X 10-" 
4.137 X 10-" 


0.2732 
0.3322 
0.3804 


3.899 X 10-'" 
1.806 X 10-"° 
4.052 X 10-"° 
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TABLE IX. The numerical data of quasi normal mode (g-mode) of the models in Table|l| {M = 0.5Mq and n(+) = 1.0) 
These data are plotted in FigW (b). 



label 


7^e (uiM) 


Im (ujM) 


Tie {u^W'/Mj 




bl 


4.514 X 10"^ 


4.849 X 10-^" 


0.5532 


5.944 X 10-"* 


b2 


3.510 X 10-2 


1.650 X 10"^ 


0.8234 


3.870 X 10-'^ 


b3 


2.737 X 10-2 


1.104 X 10"^ 


1.062 


4.287 X 10-® 


b4 


5.387 X 10-2 


7.374 X 10"'' 


0.7145 


9.781 X 10-* 


b5 


3.916 X 10-2 


4.150 X 10"^ 


1.139 


1.208 X 10-^ 


b6 


5.992 X 10-2 


5.865 X 10-* 


0.8738 


8.553 X lO-'^ 



TABLE X. The numerical data of quasi normal mode (f-mode) of the models in Table|!Il| {M — 1.2M0 and n(+) = 0.7) 
These data are plotted in Fig.0 (a). 



label 


Tie (ljM) 


Im {ujM) 






Tie {u^R^/KI^ 


Im {u^R^/M^ 


cl 


0.1333 


7.102 X 10-' 


1.126 


6.001 X 10-" 


c2 


0.1108 


5.692 X lO-"^ 


1.162 


5.967 X 10-" 


c3 


0.09569 


4.503 X 10-^ 


1.167 


5.490 X 10"" 


c4 


0.08579 


3.713 X 10-^ 


1.157 


5.008 X 10-" 


c5 


0.1503 


7.432 X 10-^ 


1.111 


5.494 X 10-" 


c6 


0.1200 


6.278 X 10-^ 


1.171 


6.125 X 10-" 


c7 


0.1020 


4.975 X 10-^ 


1.177 


5.739 X 10-" 


c8 


0.09205 


4.201 X lO-"^ 


1.167 


5.326 X 10-" 


c9 


0.1258 


6.560 X 10-^ 


1.175 


6.130 X 10-" 


clO 


0.1047 


5.163 X lO-'^ 


1.181 


5.824 X 10-" 


cll 


0.09564 


4.476 X lO--"^ 


1.172 


5.487 X 10-" 



TABLE XI. The numerical data of quasi normal mode (g-mode) of the models in Table|lJ (M = I.2M0 and n(+) — 0.7) 
These data are plotted in Fig.^ (b). 



label 


Tie {loM) 


Im (wM) 






Tie [to^Jm/M^ 


Im {uj^jR/M^ 


cl 


4.782 X 10-2 


8.763 X 10-"' 


0.4040 


7.404 X 10-"" 


c2 


3.870 X 10-2 


9.396 X 10-" 


0.4057 


9.851 X 10-"° 


c3 


2.701 X 10-2 


1.154 X 10-" 


0.3293 


1.407 X 10-"° 


c5 


5.926 X 10-2 


1.799 X 10-"° 


0.4381 


1.330 X 10-° 


c6 


4.725 X 10-2 


2.931 X 10-"° 


0.4610 


2.859 X 10-° 


c7 


3.238 X 10-2 


2.666 X 10-"" 


0.3735 


3.075 X 10-"° 


c9 


5.346 X 10-2 


5.412 X 10-"° 


0.4996 


5.058 X 10-° 


clO 


3.578 X 10-2 


3.374 X 10-"" 


0.4036 


3.805 X 10-"° 
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TABLE XII. Comparison with the results of Cowling approximation. The labels in the table corresponds to each model. 
The data in this Table is only for p(+) = 8.0 x lO" (g/cm^). 



label 




e.rror^o> (%) 




error^' ' (%) 


al 


0.2689 


-1.570 


1.304 


12.01 


bl 


0.5533 


0.01808 


1.520 


15.10 


b2 


0.8329 


1.161 


1.733 


17.94 


b3 


1.107 


4.216 


1.717 


17.41 


cl 


0.4015 


-0.6188 


1.291 


14.63 


c2 


0.4016 


-1.003 


1.355 


16.57 


c3 


0.3232 


-1.837 


1.379 


18.21 


c4 






1.388 


19.97 
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